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Abstract. Light rays from a multiply imaged quasar usu- 
ally sample different path lengths across the deflector. Ex- 
tinction in the lensing galaxy may thus lead to a differ- 
ential obscuration and reddening between the observed 

: macro-lensed QSO images. These effects naturally depend 
on the precise shape of the extinction law and on the 
redshift of the lens. By means of numerical Monte-Carlo 
simulations, using a least-squares fitting method and as- 
suming an extinction law similar to that observed in the 
Galaxy, we show how accurate photometric observations 
of multiply imaged quasars obtained in several spectral 
bands could lead to the estimate of the lens redshift, ir- 
respective of the visibility of the deflector. Observational 
requirements necessary to apply this method to real cases 
are thoroughly discussed. If extinction laws turn out to 

' be too different from galaxy to galaxy, we find out that 
more promising observations should consist in getting high 
signal-to-noise low resolution spectra of at least three dis- 
tinct images of a lensed quasar, over a spectral range 
as wide as possible, from which it should be straightfor- 
ward to extract the precise shape of the redshifted extinc- 
tion law. Very high signal-to-noise, low spectral resolu- 
tion, VLT observations of H1413+117 and MG 0414+0534 
should enable one to derive such a redshifted extinction 
law. 

Key words: gravitational lensing - Galaxies: redshifts - 
dust, extinction - Methods: numerical 



terns. Usually, the lens redshift is determined from spec- 
troscopic observations of either emission lines due to the 
lensing galaxy or of galaxy absorption lines detected in the 
spectra of the QSO images. In this work, we present an al- 
ternate method to derive the redshift of a too faint and/or 
an invisible lens. First, we describe the method which is 
based upon the fitting of photometric observations with a 
redshifted extinction model and we then show results from 
Monte Carlo simulations. The latter turn out to be very 
useful in order to assess the errors affecting the estimate of 
the lens redshift. Next, we attempt to apply this method 
to several known gravitational lens systems and outline 
all observational requirements needed in order to derive 
successful results. Finally, we propose a generalization of 
this method. Conclusions form the last section. 



2. Description of the method 

2.1. Generalities 

Figure |l| depicts a gravitational lens system. The light 
rays emitted from the quasar are deflected by the lens- 
ing galaxy and the resulting macro images are reddened 
because of extinction effects in the deflector. Due to the 
different path lengths across the galaxy, differential red- 
dening results between the macro lensed components. 



1. Introduction 

For the case of several known gravitational lens 
systems, the deflector is either not yet detected 
or it is so faint that its redshift remains un- 
known (see the Gravitational Lensing web page at 
http : //vela, astro .ulg. ac . be/gravJ.ens/). Measure- 
ment of the deflector redshift is mandatory in order to 
properly model the corresponding gravitational lens sys- 
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Observer Deflector Quasar 

Fig. 1. Schematic representation of a multiply imaged 
quasar. La — Lr> denote different path lengths across the 
deflector. 
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The flux received by the observer from component i at 
a wavelength A is: 



-Ext(A',fl v )ii 



(1) 



where F® 



is the flux that would be observed from the 
quasar in the absence of lensing, Ai the macro lens am- 
plification of component i, Ext(A',i?y) the Galaxy dust 
opacity at the wavelength A' = (z\ being the lens 



redshift) and Ry = g^L (see Cardelli et al. |1989| ) , and 
Li the path length of component i through the deflector. 

When deriving Eq. (|l|), we have of course implicitly 
assumed that the macro lensed images are not contami- 
nated by light from the lensing galaxy, that they are not 
affected by microlensing effects, that no intrinsic colour 
variations exist between the lensed components, that the 
law of extinction (first assumed to be alike that of the 
Galaxy) is similar all through the deflector and that the 
observed reddening takes place in a single lensing plane. 
A check on the applicability of most of these conditions 



can be made as suggested in section 2.3. 

The colour of component i, between the wavelengths 
Ai and A2, is then merely given by: 



m X2 



-2.5 log 



' F° e -Ext(Ai,iJv)i/ 
Al 



F° e 

A2 



-Ext(\' 2 ,R v )Li 



(2) 



Having chosen one of the macro-lensed images as a refer- 
ence component (i.e. i = ref), the relative colour of com- 
ponent j (i.e. i — j) may be expressed as: 

(mAu - m\ 3l j) ~ {m\ uTe f ~ "iA 2 ,rof) = 

C, [Ext(A' 1 , Ry) - Ext(A' 2 , R V )\ (3) 

where Cj = Lj — L m f represents the relative path length 
of component j. 

2.2. Observational constraints from the measurement of 
colours 

If we assume that the gravitational lens system consists 
of N macro lensed components and that observations are 
being carried out through M filters, we easily find that 
there are M — 1 independent colours for each of the N — 

1 components (the reference component is of course not 
included). Thus, there results a set of (M — l)(iV — 1) 
equations alike Eq. (||). 

Besides, there are 2 unknowns which do not depend 
on the numbers of filters and components: z\ and Ry. 
Moreover, for each of the N — 1 components, there is one 
additional unknown: Ci. Consequently, we have a total of 

2 + (N-l)xl = N + l unknowns. 

Since we must have more equations than unknowns, 
we easily establish the necessary condition: 



For N = 2, 3, or 4, we easily find out that M ^ 4, 3 or §, 
i.e. a minimum of 4, 3 or 3 filters is requested, respectively. 

Therefore, observational data of multiply imaged quasars 
acquired through several photometric filters are manda- 
tory. 

Note that if data obtained in M filters come from sev- 
eral different observing runs, it is then safer to only use 
those colours formed between filters from the same run 
and hence we finally get less than M — 1 useful colours. 
The above necessary condition may consequently not be 
any longer fulfilled. 

2.3. Colour- colour diagrams 

Before applying the proposed method, we must make sure 
that in a colour-colour diagram, and in accordance with 
Eq. (||), the measurements of the various lensed compo- 
nents are approximately located along a straight line. Oth- 
erwise, it may mean that one or more of the conditions 



listed in section 2.1 are not fulfilled. 

Secondly, the differential reddening for most of the 
components has to be significantly larger than the obser- 
vational errors. Otherwise, the slope of the straight line 
in the colour-colour diagram remains too uncertain. For 
example, the components of the gravitational lens sys- 
tems PG 1115+08 (Kr istian et al. |1993| ) and B 1422+231 
(Yee & Ellingson |l994j ; Yee & Bechtold |1996| ) have very 
similar colours and hence these cases are not appropriate 
for application of this method. Promising candidates are 
H1413+117 (see an example of colour-colour diagram in 
Fig. |[) and especially MG 0414 +0534 which is known to 
be very stable (Moore & Hewitt 1997) and probably not 



affected by significant microlensing effects (McLeod et al. 



1998) 
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Fig. 2. Example of a colour-colour diagram for the case 
of H1413+117, constructed from observational data listed 
in Table |^. Image A has been selected as the reference 
component. 
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Fig. 3. x 2 ma P from simulations of '5 filters' observations, 
adopting Ry — 3.1, z; = 0.5, C = 0, 0.6, 1.3, l.f for com- 
ponents A, B, C and D, respectively and a noise of 0.1 
mag. for the simulated U, B, V, R and I magnitudes. The 
symbols © and * correspond to the input and calculated 
values of the zi and Ry parameters. 
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Fig. 4. Lens redshift dispersion a versus the lens redshift 
zi for five different levels of noise and adopting Ry = 3.1, 
C = 0,0.6, 1.3, 1.1 for components A, B, C and D, respec- 
tively. The photometric observations have been simulated 
for the five U. B, V, R and I Johnson filters. 



2.4-. Fitting of the measured colours 

The left side of equation (||), which corresponds to an 
experimental relative colour, is directly derived from pho- 
tometric observations of the macro lensed quasar images. 

The right side of equation (||), and more explicitly the 
quantity "Ext(A' 1 , Ry) — Ext(\' 2 , Ry)" , is systematically 
computed for different contiguous values of Ry and zi us- 
ing the Galaxy extinction law from Cardelli et al. ( 1989| ) . 
The relative path lengths Li are fitted using a routine that 
minimizes the following \ 2 '- 



E (Experimental colour — Fitted colour) 2 
12 • ( 5 ) 



<7 



Experimental colour 



Figure |3] presents such a % 2 map obtained from a set 
of simulated data (see below). 



3. Simulations 

In order to study the errors affecting the estimate of the 
lens redshift using the fitting method described above, we 
have carried out a wide range of Monte Carlo simulations. 
When doing these, we first fix the values of the following 
parameters: the lens redshift z\ within a defined interval, 
the parameter Ry usually set to 3.1 - corresponding to the 
average galactic value -, the relative macro-amplification 
Ai = and the relative path length Li through the de- 
flector for each lensed component so that we can calculate 
precise relative magnitudes for as many different filters 
as we wish. Then we add to these magnitudes a Gaussian 
noise of a certain level and apply the above fitting method. 
We repeat these two last steps 100 times and calculate the 
dispersion a in the estimate of the lens redshift and the 
Ry parameter. 




Fig. 5. Extinction law of Cardelli et al. for Ry = 3.1. 
The redshifted (z — 1.55) positions of selected filters are 
indicated. 



Carrying out such simulations for different contiguous 
values of the theoretical lens redshift and for levels of noise 
comprised between 0.01 mag. and 0.05 mag., we obtain the 
curves illustrated in Fig. || for which we used the set of 
the five U, B, V, R and I Johnson filters. We see that the 
dispersion a affecting the estimate of z; is about 0.1 and 
even less at higher redshifts. The visible "oscillations" in 
Fig. U correspond to the positions of the wavelengths of 
the different selected filters with respect to the redshifted 
2175A absorption bump (see Cardelli et al. 1989). As an 
example, Fig. |5| represents the extinction law of Cardelli 
et al. ( |1989| ) as a function of the reciprocal wavelength j- 
for Ry = 3.1. The redshifted (z = 1.55) positions of the 
U, B, V, R, I, J, H, K and L filters are indicated. 
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Fig. 6. Histograms of the repartition of Ry for two dif- 
ferent levels of noise and adopting Ry = 3.1, zi from to 
2.65 and C = 0,0.6, 1.3, 1.1 for components A, B, C and 
D, respectively (see Fig. |Z]). 

Furthermore, we have illustrated in Fig. ^| dispersions 
affecting the derived values of Ry for the cases of 5 filters 
and two different levels of noise: 0.01 mag. and 0.05 mag. 
We see that, in this latter case, the dispersion gets rather 
large. 

Next, Fig. ^ shows on a logarithmic scale differences 
in the lens redshift dispersion when using the U, B, V, R 
and I filters or the infrared I, J, H, K and L filters. Since 
the method is sensitive to the extinction gradient between 
two filter wavelengths, the much higher dispersion derived 
for the infrared filters is mainly accounted for by the shal- 
lower shape of the extinction law in the infrared range. 
Moreover, the lower extinction in the infrared range also 
affects the dispersion in the estimate of the lens redshift 
because the experimental relative colours are smaller and 
hence the effects of noise get more pronounced. 

Finally, we have also used other extinction laws such 
those proposed by Pei (Pci 1992) for the Milky Way 
("MW"), the Large Magellanic Cloud ("LMC") and the 
Small Magellanic Cloud ("SMC" ). From the same set of 
simulated data as in Fig. y, Fig. |] shows the resulting x 2 
curves. Comparable minima and lens redshift solutions are 
usually found for the MW and LMC laws. However, high 
quality photometric observations are mandatory in order 
to efficiently select the most likely lens redshift among ac- 
ceptable solutions for the different extinction laws. Note 
that because of the noise in the simulated photometric 
data, the minimum for the MW curve in Fig. || is found 
at zi — 0.48, instead of the theoretical value zi = 0.5. 

4. Applications to some real gravitational lenses 

4.1. 2237+0305 

Figure || represents the x 2 map for the Einstein Cross de- 
rived from the published sets of photometric observations 



Fig. 7. Lens redshift dispersion a versus the lens redshift 
zi for two different sets of filters and adopting Ry =3.1, 
C = 0, 0.6, 1.3, 1.1 for components A, B, C and D, respec- 
tively. A noise of 0.01 mag. has been used for the simulated 
photometric observations. 




Fig. 8. x- 

tions adopting three different extinction laws and Ry = 
3.1, zi = 0.5, C = 0,0.6,1.3,1.1 for components A, B, C 
and D, respectively and a noise of 0.1 mag. for the simu- 
lated magnitudes (see text and Fig. |^) . 

listed in Table § The minimum is found for zi ~ 0.13. 
However, the redshift of this gravitational lens is known 



to be 0.04 (Huchra et al. 1985). The small discrepancy 
between these two redshift values is very likely due to the 
fact that the data were obtained during two different runs 
when photometric quality was not good enough (ct ^ 0.02 
mag. for approximately half the reported values in Table 
Q). Moreover, this gravitational lens system is very likely 
affected by microlensing effects on a time scale of several 



years (Irwin et al. 1989; Corrigan et al. 1991; 0stensen et 
al. 1996j ). The contours in Fig. [| correspond to Ax 2 = 1 
and delineate, at a confidence level of 68%, acceptable 
ranges for the values of Z[ and Ry. Very high signal-to- 
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Table 1. Photometry of 2237+0305 

Components 



Date 
25/09/87 
25/09/87 
25/09/87 
10-11/10/95 
10-11/10/95 
10-11/10/95 
10-11/10/95 



Filter 



i 

B 
V 
R 
I 



A 

17.77±0.03 
17.62±0.02 
17.25±0.02 
17.626±0.009 
17.278±0.004 
17.093±0.004 
16.966±0.006 



B 

17.98±0.03 
17.77±0.03 
17.38±0.02 
17.741±0.008 
17.428±0.009 
17.274±0.004 
17.192±0.007 



C 

18.46±0.05 
18.06±0.04 
17.61±0.03 
18.881±0.020 
18.389±0.009 
18.109±0.008 
17.910±0.015 



D 

18.69±0.06 
18.40±0.06 
17.98±0.04 
19.111±0.033 
18.734±0.022 
18.441±0.010 
18.206±0.020 



Reference 3 
1 
1 
1 
2 
2 
2 
2 



(1) Yee 1988; (2) Burud et al. 1998 (CLEAN algorithm 




Fig. 9. x 2 map for 2237+0305 derived from the published Fig. 10. \ 2 map for H1413+117 derived from the pub- 
sets of photometric observations listed in Table fjj The x 2 lished sets of photometric observations listed in Table |[ 
minimum corresponds to z\ = 0.13 and Ry = 2.8. The x 2 minimum corresponds to zi = 1.15 and R v ~ 9. 



noise photometric observations of the Einstein Cross, us- 
ing several filters during a single run and reliable tools to 
derive the photometry of the multiple QSO images, are 
badly needed to confirm that the different experimental 
relative colours reported between the lensed components 
arc effectively, and solely, due to extinction effects in the 
deflector. 



4.2. H1413+117 

Figure [l(] shows the x 2 map for the Cloverleaf. For this 
case, we used photometric data obtained through 7 filters 
but from two distinct observing runs (see Table |J) . We see 
that a minimum is found for zi ~ 1.15 and Ry — 9. This 
latter value looks quite unexpected and we also conclude 
here that very high signal-to-noise photometric data of 
the Cloverleaf, using several filters during a single run, 
are badly needed to confirm that the reddening observed 
between the macro-lensed QSO images is essentially due 
to extinction effects. For this gravitational lens system, 
numerous lens redshifts from identified absorptio n line s 

|l996p , 



have been proposed: 0.31 and 0.75 (Afanas'ev et 



al. |1984| ; Drew & Boksenberg |l984j Turnshek et al. [1988|) 
but the real lens redshift still remains unknown. 



4.3. MG 0414 + 0534 

Figure |ll| shows the x 2 map for MG 0414+0534 derived 
from the published sets of photometric observations listed 
in Table ||. Because we only have 2 colours, including the 
near infrared colour H— K (see Table|3|), the estimate of the 
lens redshift by this method remains very uncertain (see 
the very widespread contour in Fig. p+| ). Due to the known 
flux stabil ity of MG 0414+0534 in the radio (Moore & He- 
witt 1997 ), this gravitational lens system constitutes one 



of the most promising targets for an estimate of the lens 
redshift from the observed reddening of its lensed compo- 
nents. Microlensing eff ects m ust be insignificant for this 
object (McLeod et al. 1998 ). High quality and simulta- 



neous photometric data from the visible up to the near 
infra-red are very much needed for this interesting sys- 
tem. 



1.44 (Magain et al. |j988| ), 1.66, 1.87 and 2.07 (Hazard et 
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Table 2. Photometry of H1413+117 



Components 



Date 


Filter 


A 


B 


C 


D 


Reference 3, 


02/08/94 


B 


-1.06±0.02 


-0.82±0.02 


-0.88±0.02 


-0.70±0.02 


1 


02/08/94 


V 


-0.41±0.01 


-0.23±0.01 


-0.20±0.01 


-0.06±0.01 


1 


02/08/94 


R 


-0.08±0.01 


0.10+0.01 


0.18±0.01 


0.28±0.01 


1 


02/08/94 


I 


0.19±0.01 


0.36±0.01 


0.48±0.01 


0.62±0.01 


1 


22/12/94 


F336W 


19.03±0.03 


19.37±0.02 


19.01±0.01 


19.20±0.03 


2 


22/12/94 


F702W 


18.40±0.01 


18.55±0.02 


18.68±0.01 


18.75±0.01 


2 


22/12/94 


F814W 


18.66±0.02 


18.76±0.02 


18.96±0.02 


19.06±0.02 


2 



(1) 0stensen et al. 1997; (2) Turnshek et al. 1997 



Table 3. Photometry of MG 0414+0534 



Date 


Filter 




Components 




Reference 3 


A 


B 


C 


D 


10/03/93 


H 


15.31±0.05 


15.51±0.05 


16.28±0.01 


17.16±0.02 


1 


10/03/93 


K 


14.34±0.05 


14.5640.05 


15.3640.01 


16.24±0.02 


1 


08/11/94 


F675W 


22.760±0.008 


23.75640.016 


23.48840.012 


24.258±0.022 


2 


08/11/94 


F814W 


20.59540.003 


21.40740.004 


21.36340.004 


22.212±0.007 


2 



(1) McLeod et al. |1998| ; (2) Falco et al. |1997| . 




Fig. 11. x 2 map for MG 0414+0534 derived from the pub- 
lished sets of photometric observations listed in Table 
The x 2 minimum corresponds to zi = 1.15 and Rv = 1.8. 

5. Generalization 

From Eq. ([!]) and considering the j component relative to 
the re/(-erence) component, we may also write: 



r \ _ J}±_ „-Ext(\' ,R v )[L j -L IQ{ ] _ a -Ext(A', R v )Cj 

F rcf ~ A f ~^o e 



(6) 



where Aj = . From this equation, we straightforwardly 
derive: 



m\j - mA,ref 



2.5 
ln(10) 



Cj [Ext(A', Ry)] - 2.51og(A, 



(7) 



We see that in fact, one can also directly work with the 
relative magnitudes m\ j — ^x,rcf 

but that the relative 
macro lens amplifications Aj then appear as additional 
unknowns which also need to be fitted. Following a rea- 
soning similar to that adopted in section [2.2| , the neces- 
sary condition is also found to be: M(N — 1) ^ 2N, where 
M and N have the same meaning. Before applying this 
generalized method to real photometric observations, one 
should make sure that intrinsic or extrinsic (cf. induced by 
micro-lensing) light variations are not affecting the indi- 
vidual lensed components. These conditions are of course 
more difficult to be checked than those imposed on the 
observed colours. 

Finally, it is also straightforward to establish from Eq. 
(0) that the redshifted extinction law Ext (A') may be di- 
rectly retrieved as a function of A' from the observed spec- 
tra of at least three multiple lensed QSO images. Very 
high signal-to-noise, low spectral resolution, VLT obser- 
vations of gravitational lens systems like H1413+117 and 
MG 0414+0534 should enable one to derive such a red- 
shifted extinction law and compare it with the known 
MW, LMC and SMC ones. 



6. Conclusions 

Monte Carlo simulations show that we need very accu- 
rate photometric data (typically a ~ 0.01 mag.), in sev- 
eral (^4 — 5) filters, in order to estimate the lens red- 
shift from the differential reddening observed between the 
lensed images of a multiply imaged quasar. Furthermore, 
the macro-lensed images must not be affected by light con- 
tamination from the lensing galaxy, by microlcnsing effects 
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or by intrinsic colour variations reflected with a time de- 
pendence between the lensed components. Unfortunately, 
all these conditions are rarely fulfilled at the same time 
and existing photometric data are often not good enough 
or not sufficiently reliable. Consequently, accurate pho- 
tometric observations of gravitational lenses obtained in 
several filters during the same run are badly needed. Fi- 
nally, we have seen that the redshifted extinction law of 
the lensing galaxy could be directly retrieved by means of 
a similar method provided that we may record the indi- 
vidual spectra of at least 3 lensed components, over a wide 
wavelength range, assuming that no light contamination, 
no microlensing effects, . . . affect the individual lensed 
components. Very good such candidates to be observed 
with the VLT are H1413+117 and MG 0414+0534. 
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Note added in proof: On the basis of spectroscopic ob- 
servations obtained for MG 0414+0534 w ith the KECK 
telescope (see preprint astro-ph / 9809063] submitted on 5 
September 1998), J.L. Tonry and C.S. Kochanek report 
the redshift determination (zi ~ 0.96) of the very faint 
gravitational lens. Based upon preliminary photometric 
observations of this system, our own redshift estimate is 
found to be in reasonable agreement with this value (see 
Fig. IB). 
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